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TREATMENT OF THE CONTROL MECHANISMS 0? LIGHT AIRPLANES 

IN THE FLUTTER CLEARANCE PROCESS 

Elnar J. Breitbach* 

Langley Research Center 


SUMMARY 

Recently, it has become more and more evident that many difficulties 
encountered in the course of aircraft flutter analyses can be traced to strong 
localized nonlinearities in the control mechanisms. To cope with these prob- 
lems, more reliable mathematical models paying special attention to control 
system nonlinearities may be established by means of modified ground vibration 
test procedures in combination with suitably adapted modal synthesis approaches. 
Three different concepts are presented in detail. 


INTRODUCTION 

At first glance the flutter clearance of soaring and light airplanes does 
not seem to raise any serious problems which cannot be solved by means of 
today’s aeroelastic tools. This is true even for the determination of the 
unsteady aerodynamic loads as long as cases with large aspect ratios at compa- 
rably low speeds are considered. The elastodynamical characteristics can be 
determined by using common experimental or analytical methods if structural 
linearity can be assumed to be a proper approximation. However, as experience 
has shown, the control mechanisms of light airplanes 1 are generally nonlinear 
to such a large extent that setting up a dependable mathematical model requires 
special attention, including modifications to standard linearized procedures. 

In the first part of this paper some of the most frequently occurring types 
of control-system nonlinearities are described. To get an idea of the influence 
of some typical nonlinearities on the aeroelastic stability the results of wind 
tunnel flutter tests on a nonlinear wing aileron model are presented. After 
that, it is shown in detail how the aeroelastic equations of light airplanes 
with localized nonlinearities may be formulated by using various suitably modi- 
fied ground vibration test (GVT) procedures all based on the well-known modal 
synthesis approach. The shortcomings as well as the usefulness of the different 
concepts are discussed. 


♦NRC-NASA Senior Resident Research Associate. 

1 Light airplanes as sued in this paper include both powered and unpowered 
vehicles where the power to the flight control system is supplied by the pilot 
without electrical or hydraulic boost through a system of cables, pulleys, push- 
rods, bellcranks, or other mechanical linkages. 
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SYMBOLS 


a»b 

A,B,C 


hinge axis coordinates of oontrol surfaces and tabs, respectively 

stiffness Matrices, respectively, defined in terns 
of geometrical displacements 


Afc,AB,Ac “trices of mass, damping, and stiffness changes, respectively, 
defined in terns of geometrical displacements 

Ce 
e»f 


equivaienti inear stiffness of a nonlinear force deflection diagram, 
defined in equations (1) and (2) 

center-of-gravity coordinates of control surfaces and tabs, 
respectively 

force or moment acting on a control surface or tab 
column matrix of constraint functions gj 

bending deflection of the quarter-chord line of lifting surface 

mass moments of inertia per span unit of control surface and tab, 
respectively, referred to the center of gravity 

I Rxf 1 Ry»lR Z «ass moments of inertia of control surface referred to the main 
axes c* inertia 


F 

g 

h 

Iv.Itv 


l 

®R 

■v»®tv 

M,D,K 


span width coordinate 
control- surf ace mass 

mass per unit span of control surface and tab, respectively 
generalized mass, damping, and stiffness matrices, respectively 


4M ' iD ’“ lnto * ccou " t **•■• - •**«— 


p 

g r p 

Q 

t. 

T 

u 

0 


column matrix of external forces 
column matrices of generalized coordinates 
column matrix of generalized forces 
time 

inertia energy 

column matrix of geometrical deflections 
stiffness energy 


v flight speed 

W damping energy 

X,y transformation matrices, defined in equations (53) and (55) 

a rotation about the quarter-chord line of lifting surface 

0 control-surface rotation about the hinge line 

Y tab rotation about tab hinge line 

n rotation of a control surface referred to its center of gravity 

0 damping loss angle 

X column matrix of Lagrange's multipliers H 

A diagonal matrix of the square values of the normal circular 

frequencies 

^ modal matr ices 
M circular frequency 

I unity matrix 

0 zero matrix 


j « imaginary unit 

Subscripts t 


A,B,C,R,v,t substructures indices 
j constraint index, l • 1, 2, ...» o 

t normal mode index 

a , e indices referring to o constraints and e independent coordinates 

NL index referring to nonlinear properties 

L index referring to linear properties 


Superscripts: 

<j transposed matrix 

A,B indices referring to substructures A and B 

real, imaginary part 


i 


4 
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GENERAL REMARKS 




Sources of Control-System Nonlinearities 

Aeroelastic investigations are usually carried out on the basis of simpli- '• 

fled linearized mathematical models. In many cases this approach has been ade- 
quate to ensure sufficient flutter safety margins for light airplanes. However, 
in the last few years, it has become evident that disregarding nonlinear phenom- 
ena can lead to hazardously misleading results. For example, it is shown in 
reference 1 that so-called concentrated or localized nonlinearities in control 
systems have a significant effect on the flutter behavior. Nonlinearities of 
this kind may be produced by such things as 

(1) Backlash in the joints and linkage elements 

(2) Solid friction in control-cable and pushrod ducts as well as in the 

hinge bearings 

(3) Kinematic limits :ic*r. of the .ontrol-surface stroke 

(4) Application of special spring tab systems provided for pilot handling 

relief 

The most critical parts of a control mechanism where localized nonlinearities 
may arise are shown schematically in figure 1 . 

An aeroelastic investigation may become even more complicated if it is 
necessary to account for items such as the following: 

(1) Preload changes due to maneuver loads and specially trimmed flight 

attitudes 

(2) Changes in friction and backlash over an airplane's lifetime 

(3) Additional mass, stiffness, and damping forces randomly activated by 

the pilot 

Coping with all these difficulties requires special measures throughout the 
flutter clearance process. First, the ground vibration teft (GVT) used to 
determine the elas todynamical coefficients of the flutter equations has to be 
modified so that a consistent and superpositionable set of orthogonal, or 
well-defined nonorthogonal, normal modes can be measured. 

In reference 2 a proposed experimental approach employs a high frequency 
auxiliary excitation superimposed upon the much lower sinusoidal excitation to 
be tuned to the several normal frequencies. Thus, "slip-stick" effects and 
related nonlinearities in the control mechanisms can be minimized. The method 
requires additional test and control devices capable of sxciting all controls 
simultaneously. 

Of course, the simplest solution appears to be to build control -surface 
mechanisms without either friction or backlash. However, aside from a consider- 
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able increase in manufacturing costs, there is no guarantee that such an ideal 
condition could be kept unchanged for the lifetime of an airplane. Moreover, 
a fr ictionless control system is not necessarily equivalent to better handling 
qualities, because friction helps give the pilot the "feel" of flying the 
airplane. 

From an experimentalist's standpoint, there are some simpler, but effective, 
methods using special modal coupling and modal superposition approaches. A 
detailed presentation of some of these methods is given in the subsequent sec- 
tions of this paper, they will be referred to as Concepts I, II, and III. 


Illustrative Examples of Control-System Nonlinearities 

To get a realistic impression of control-mechanism nonlinear it ies, the 
force deflection diagrams F(fJ) of the rudder and aileron system (anti symmet- 
rical and synmetrical case) of a soaring airplane (ASW-15, A. Schleicher, 
Poppenhausen, W. Germany) are shown in figures 2(a), 3(a), and 4(a). Using the 
principle of the energetic equivalence (refs. 1 and 3) the stiffness and damp- 
ing properties of a nonlinear force deflection diagram can be approximated by 
the so-called equivalent complex stiffness: 


C e (B) - Ce(B) + j Ce(B) 


(j - (i) 


• « 

The coefficients C e (8) and C e (8) , representing stiffness and damping, respec- 
tively, can be calculated from 


, 1 f 2 " 

C„(B) ■ — \ F(B cos 4>, -Bio sin <|>) cos <J> d<p 

*8 \=Q 

n 1 

C-(3) ■ — \ F(B cos <J>, -Bw sin <t>) sin $ d<|> 


( 2 ) 


where the circular frequency w * 2nf (where f is frequency in hertz) and 
the integration variable <)> ■ wt. Damping can also be expressed by the loss 
angle 


N 


C e (8) 



( 3 ) 
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The functions C e (B) and C e (@) corresponding to the force deflection diagrams 
of figures 2(a), 3(a), and 4(a) are plotted in 2(b), 3(b), and 4(b), respec- 
tively. Figure 3(b) shows that the antisyanetric aileron hinge stiffness in the 
range of the normal aileron stroke varies between 390 N-m and 44 N-m. Because 
of the stiffness variation, the normal frequency of the antisynaetrical aileron 
vibration (wing assumed to be fixed) varies over a wide range, between 2.4 Hz 
and 7.4 Hz. At least two other antisyanetric normal nodes lie in this frequency 
range and are consequently characterized by highly amplitude— dependent portions 
of aileron vibrations. Similar effects can also be observed for the symmetric 
aileron aode and for the rudder vibration. 

The effects of strong nonlinearities on the flutter behavior have been dem- 
onstrated in some wind-tunnel tests carried out on a nonlinear wing-aileron model 
in the low-speed wind tunnel of DFVLR Gottingen. The nonlinear flutter bound- 
aries for a backlash-type and for a spring-tab-type aileron hinge stiffness are 
shown in figure 5. Unlike the flutter boundaries of linear systems, both curves 
are characterized by a considerable dependence of the critical flutter speed on 
the aileron amplitude. Thus, the flutter boundary of the spring- tab- type system 
varies between V * 12.5 %'s and V * 24 m/s. The backlash- type system shows 
a flutter boundary variation between V ■ 13.5 m/s and V = 20 m/s. More 
detailed information, especially about the geometric and elastodynamic data of 
the wing-aileron model, is presented in reference 1. 


MATHEMATICAL MODELING USING MODAL SYNTHESIS CONCEPTS 

As mentioned previously, the determination of the elastodynamic character- 
istics by means of GVT can be affected severely by localized nonlinearities in 
the control mechanisms. It will be shown in the following discussion that the 
uncertainties resulting from these nonlinear effects can be avoided by applying 
experimental-analytical concepts based on the well-known modal synthesis 
approach. 

Each of these concepts can be used to set up the aeroelastic equations of 
the actual airplane including all control-mechanism nonlinearities. The non- 
linear force deflection diagrams of the different controls can be determined 
by static or dynamic tests. 

Three different concepts will be presented. They may be briefly described 
as follows: 

Concept I: Measurement of a set of orthogonal normal modes with the con- 

trol surfaces rigidly clamped; separate determination of the control- 
surface normal modes with the rest of the airplane rigidly fixed. 

Concept II: GVT on a configuration artificially linearized by replacing 

the nonlinear control- mechanism elements by linear and lightly damped 
dunmy devices; thus, a set of orthogonal normal modes for the entire 
system is available. 


♦ 
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Concept Ills Measurement of a set of orthogonal normal modes with the 
control surfaces removed; separate determination of the normal modes of 
the control surfaces in uncoupled condition. 


Concept I 

The governing equations of motion of an aeroelastic system, formulated in 
terms of physical coordinates, can be written in matrix notation as follows; 


Au + Bu + Cu 


(4) 


where 

A mass matrix 

B damping matrix 

C stiffness matrix 

u column matrix of the physical displacements; u and ii are first 

and second derivatives with respect to time t 

P column matrix of external forces, for instance, unsteady aerodynamic 

forces 

It is obvious that parts of the matrices B and C are nonlinear because of 
the localized nonlinearities of the controls. 


Controls without tabs.- if the GVT is carried out with the controls rigidly 
c amped to the adjacent structure, a set of n largely linear normal modes 
^Ar d® measured and combined in the modal matrix 


*A * [*A1» *A2* • • •» *Ar» • • 



(5) 


The modes satisfy the orthogonality condition 
♦a^a “ m A 'l 

*a t »a - A A M A - K A J (6) 

where 

diagonal matrix of the generalized masses 
K \ diagonal matrix of the generalized stiffnesses 
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a A diagonal matrix of the square values of the circular normal 

frequencies 10 ^ 

The generalized damping matrix d a (not necessarily diagonal) is defined by 
°A « V^A 


iS«eS S ^^r^ e ^ r °\ S “ faCeS are ri9id in the ^ re< 3 uenc y range of 

f additlonal control-surface rotation modes with the 

ln 8 ructure at rest can be determined and combined in the modal 


[*B1r *B2* ...» «>BVr . . OgJ 


SLS^rdinate^br 0 ' 5 ° f ^ CO,nPlete Structure are rel ^ed to the gener- 


U ■ 4>q 


where the column matrix of the generalized coordinates 
* “ &A T ' < 3b T ] T 


q r and is 


* ■ ♦»] ni> 

The basic idea of this modal superposition is outlined in fiaure 6 « 5 nh«n <-..«■ 
ing equation (9) into equation (4) and premultiplying it by 4 >T yields 


Mq + Dq + Kq « Q 


where 
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\ ( 13 ) 

°1 

4, 

The matrices M^, K^, and measured in a GVT are defined in the equa- 

tions (6) and (7). The diagonal matrices Mg, Dg, and K B contain the gener- 
alized masses, damping values, and stiffnesses of the control- surf ace rotation 
modes. In the case of nonlinear hinge stiffness and damping, the matrix elements 
of K b and D B are 

« 2 c ev(f$) 2 

Kjjv = C e >j(B) B\jq Dbv = — (14) 


• " 

where C €V (&) and C eV (B) can be determined from equation (2). The term B V a 
denotes the control rotation in the action line of the control actuator force. 
The matrix can be determined by calculation or measurement caking into 

account not only the con trol-sur face mass but also the moving mass of such 
attached hardware as pushrods, cables, and control stick. The elements of the 
coupling matrix 

m AB = m BA T = *A T A*b ^ ^ 


Ma 

Mab 

Mba 

Mb 

"Da 

0 

0 

Db 


Q » 4> T P 


K = 


K A 

|0 


can be found by integration over surfaces S Bv of the controls 

M AB,rv * ^ + pv + ©\l n v( a v + dfl- (16) 

Sbv 

where the following terms correspond to the vth control with tab locked to the 
control 


m v mass of the control surface per unit span 

I sj mass moment of inertia per unit span referred to the center of gravity 

e v distance between center of gravity and hinge axis (see fig. 7) 

a v distance between hinge axis and the quarter-chord point (see fig. 7) 

All these data as well as the amplitudes h c , ct r , and B v (see fig. 7) are 
functions of the span coordinate fc. In case of an ideal locking of the con- 
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trols, neither hinge stiffness forces nor hinge damping forces are generated in 
the normal modes $Ar* *»ence. 


k AB = k ba t - 
d AB = °ba t - 


(’7) 


Extension to controls with tabs .- The above procedure can easily be 
extended to systems with controls and tabs (spring tabs, trim tabs, or geared 
tabs) by introducing the tab movement as a separate degree of freedom. For 
this special case the main GVT configuration is characterized by controls 
locked to the adjacent airplane structure and tabs locked to the controls. 

This leads to the same set of normal modes $ A r as defined in equation (5). 
Furthermore, the degrees of freedom of the controls are separately determined 
with the main structure at rest and with tabs locked to the controls. The 
resulting normal modes are identical to the ones defined by equation (8) . 
Finally, in o third step the tab inodes 4>cv ara determined with both the main 
structure and the controls at rest. This concept is schematically illustrated 
in figure 8. In accordance with this, u can be expressed as a series expan- 
sion of the normal mode sets $b* ar *d 


u - ($ A , <J> B » * C ] <1 (18) 

where 

q = (q B T . q B T . qc T ] T < 19 > 


Replacement of u in equation (4) by equation (18) and premultiplication by 
4 > t leads to an equation similar to equation (12). Because of the additional 
tab degrees of freedom the matrices M, D, K, and Q have the extended form 



Ma 

Mab 

m ac 




] 

M = 

m ba 

h b 

m bc 

Q = 4>Tp 




_^CA 

“cb 

Me 






da 

0 

0 


k a 

0 

o" 


D » 

0 

d b 

0 

K * 

0 

k b 

0 



0 

0 

°c_ 


0 

0 

K c 

> 


(20) 
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The matrices M^, Mg, * Mgg r , K*, Kg, D*, and Dg are identical to 

the matrices defined in equations (6), (7) r (14), (15), and (16). The matrices 
an< ^ D c can be determined in the same way as Kg and Og by measuring 
the nonlinear force deflection diagrams of the tabs and using equation (2) to 
calculate 


*Cv - C eV (Y) Yva 


D cv “ Cy a (Y) Yva 


( 21 ) 


The term Yva denotes the tab rotation in the line where the force acting on the 
tab is applied. Ihe matrix Me can be determined by test or calculation. The 
ele&ents of the coupling matrix 


M AC ■ Mca t " 


( 22 ) 


can be found by integration over the tab surface 

m ac *rv * ^ { f v m tv h rYv + [jtv * f\>®tv(bv + f °*rY^ dl^ (23) 

S CV 

where the following terms correspond to the vth tab (part of the vth control) 
mtv mass of the tab per unit span 

I t v mass moment of inertia per unit span referred to the center of gravity 

fv distance between the tab hinge axis and the tab center of gravity (see 

fig. 7) 

by distance between the tab hinge axis and the control hinge axis (fee 

fig. 7) 

The quantities ®tv» f\>» ®nd ** well as h c , a r , and Yv (see 

fig. 7) are functions of the tab span coordinate £ t . The elements of the 
coupling matrix 


m BC * McB 1, ■ (24) 


between the control surfaces and the appertaining tabs can be calculated by 
integration over the tab surface Scv 

**BC,vv ■ j [Itv + ®t*v(*\> + f v >] 6 v Yv <&t (25) 

8 Cv 
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the n °" !Tal u T 0de3 *Ar can be measured with ideally locked control 
and tab binges, neither hinge damping forces nor hinge stiffness forces are 
generated in * Ar . This leads to 


o 

III 

"s 

H 

3 

Kbc * Kcb T = 

o 

III 

> 

H 

3 

d bc * d cb t = 


Concept II 


(26) 


*' ^ th * re P lacen>ent <* the control non- 
by artlflclal llnear stiffnesses results in a modified linearized 
test configuration represented in matrix notation by 


A[,u + B^u + C L u * p 


(27) 


which is formulated in terms of physical displacements. The governing dynamic 

equation* (4) bT n °" linear systen can written in the sam^ tom as 

equation (4) by subdividing the matrices A, B, and c as follows: 


A - A l - Aa l + Aa nl > 
B - B l - Ab l + ABnl J> 
C - C L - Ac l + 


tetnln th^tlfJ&i 7T 6 diffeienCe in the mass distribution 

J “ " ti “ cial linaar sy «t«n and the real nonlinear system, Ab l and 

•Sm-Sf* A? ^ T “ nd 8tiffn « 88 Properties of the artificial linear 

elements; ABn L and Ac^, describe the damping and stiffness properties of 

victim n ° nlin ? ar elements * Development of the arbitrary Secernent 
.. t0 . . ln a 8er * es expansion of the measured normal modes * Lr of the 
linearized system yields Lr 


u - * L q 


(29) 


mserting this modal transformation into equation (4) , premultiplyina bv 

,2 " r * ,uU * *" 9.r,.;.U^ ,iu.u^. o( 2.L, 

SaSilT ( ’ 2 ’’ but “ lth th * a ** pi " 5 ' •- 
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M * M L - AM l + AMnl 
D - D l - dD L + tf> NL ) 


* ■ *L ” ^*L + ^Ki 


ML. 


( 30 ) 


Furthermore? ^ Dl " *"* Kl “* “ ea8ured in a GVT on the linearized system. 


( 31 ) 


AMhl - a*L - *L T (AhNL - AA l ) 4> l *' 

adnl - “ AB l )* l > 

m nl - AK l » * L T (AC NL - tc h )* h ^ 

2 rr^TSSi. rot the '' lh 


^VL * [Bvl. ®\»2» • • •* 8\>r* • • •» $vnj 
and ABjjl - Ab l and AC NL - dc L degenerate to the 1 x i matrices 


( 32 ) 


C e (B) "I 
ABjjl - AB L B L 


o> 


^HL * AC l ■ C e (&) - C L 


( 33 > 


where the nonlinear stiffness and damping value- r' /Ay 

determined again by applying equation (21 to thl «- e and C e (&) can be 
j.*,..., 1 iuy equation U) to the measured nonlinear force 

S Jl a5r 

>*. noilinwtXJSTil.S'S'thf JSI« U r ' £etrln « "*»»•« ««• 

ficl.l linearization, to ZTio4 .m. e ° r " U “' i " 5 f '°* «» 
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Concept III 


«f aeroeiaatic equations of an airplane can also be established by means 

86 °f n ? rma m ° des “« a8ur « d in a «VT with controls removed and riqid- 
body and seme elastic normal modes of the several controls (see fig. 10) deter- 

^* XPer 1 1 L entally ° r by fairly 8i " ple calculations. The equations of motion 
of the coupled system can be set up by means of Lagrange's equations 


where 


3 /3t \ 3U 3w ^ 

»Ur/ + 3q c + 3^ “ 2 H 3^ 

Jo * 1 


(34) 


• • a# fl£ + 1 # 


* ng # 


n i 9 • • . ) 


2T = Mq 2 
2U * Kq 2 + AK q 2 ) 


(35) 


2W = Dq 2 + AD q 2 


The matrices AK and AD in equation (35) take into account the elastic coupl- 
ing between control surfaces and main structure by means of the real hinge stiff- 
ness and hinge damping elements. The term on the right side of equation (34) 
is formulated in terms of Lagrange's undetermined multipliers u which corre- 
spond to a number of o constraint conditions 1 


(<Jl * <12* • • •* <J r » • • •* * 0 01 


2 * • • • f o) (36) 


They express compatibility in those coupling points, where the controls can be 
assumed to be rigidly fixed to the main structure. Application of equation (34) 
to equations (35) and (36) yields H ' ’ 


Mq + (D + AD) q + (K + AK) q - - 0 


(37) 


where the elements of the r « l matrix f T are 
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(38) 


„ T 391 

*ir T - r“ 

r 


f Urther j deriVati ° n to the coaling of only two systems, A and 

■ass itiffoes."!^ C °" tr ° 1 ?“ face) results in the following generalized 
■ass, stiffness, and damping Matrices of the uncoupled system: 


■-[: 

-E 


*a 

0 


0 

*B 

o “ 

*8 


where the submatrices are 


°a 

o 




0 

Db 


(39a) 


Mi - 4iTAi*i Kj - ♦iTCiti Di . * i T Bi * i (i - A , B) 


(39b) 


subsystems*^ A Jnd ?!!*♦. Ci d * 8cribe m f ss ' dam P in 9' and stiffness of the 
suosystens A and B in terms of geometrical coordinates; *j is the modal 

1 ‘ Th * ele,Dents ot th « diagonal matrices «i and Ki and 

Sinld i ? e8 .w i# WhiCh ar * n0t nece8Sar ily diagonal, can be deter - 

ined by GVT or, as in the case of the controls, by calculation, also. 

can ’ *“* « — » 


dK • *AB T CAB*AB 


ad - #aB Tb AB*AB 


(40) 


Whe : i the main structure and control surface are coupled by one single complex 
hinge stiffness C e in the action line of the control foSrce, wj obtain 



A A 

A 1 


-i 

♦ab * 

a al * a a2» • • 

| • 
* 

P 
1 00 
l 3 

0 



0 

J a a,n+l » 

~B 

a a,n+2' • 

• •' a tn B 



t 
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C AB * c e 


(42) 




“l 

-f 

c m 

*1 

-T 



b ab * 





0) 



-1 

1_ 


-1 

— 


A B 

Die angles of rotation ot ac and ot ar are defined in figure 9. Por the special 
case of coupling two systems A and B the compatibility condition for o 
physical degrees of freedom can be expressed by the constraints 


gj^ x — ujk *0 (A * 1 » 2, . . . , o) (43) 

A Q 

If uj and U£ are expressed in a series of the normal inodes of the systems 
A and B, then 


n A 

9A * ^ u Ar9r 
r«l 


I 


B 

u Ar9r 


0 


(44) 


V 


or in matrix notation 


g * fq ■ 0 < 45 > 


The aeroelastic equations of motion are defined now by the n* + n^ - a gener- 
alized coordinates. Due to the o constraints there remains a number of 
e ■ m - a independent generalized coordinates in terms of which the aero- 
elastic equations have to be formulated. To do this, the term in equa- 
tion (37) has to be rearranged rowwise so that 





(46) 


where Yg is a nonsingular o * o matrix* The matrices Mr Kr Dr AK# and 
AD with respect to both their columns and rows and the column matrix q have 
to be rearranged in the same sense. The rearranged equations can be written as 

Mq + (6 + &5)q + (K + $K)q - Y T \ (47) 
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where 


5 - [p T . <te T ] T 


(4<) 


The new structure of the matrices M, 6, K, 
following equation using M as an example 


Ao, and AK is shown in the 


M ee M ea 


*oe Moo 


(49) 


Thus, X can be determined as follows 


X “ ^(Moe» Moojq + ((p ae , D ao ] + (AD oe , AD a( £|)q 

+ ([ K ae» K oo] + • ^oo])^} 


(50) 


Fra# equations (45) and (46) it follows that 


« -to^cP 


(51) 


Inserting equation (50) into the first e rows of equation (47) and taking into 
account equation (51) results in the following equation 


(**€£ “ MggX - xTMoe + X% 0 X)p 

+ { D ee + AD ee * < D eo + AD e o) x “ xT (%e + Aty, e ) + x T (D acJ + Al>t, a )x}p 
+ {k E6 ♦ AK ee - (K eo + AK ea )X - xTa^e + AK oe ) + xT^ + AK ao )x}p - O 

(52) 


where 


X - y 5 l«p c 


( 53 ) 


It can easily be shown that equation (52) can be transformed to the more con- 
venient equation 


Y*ft*p + Y^D + A5)Yp + Y t (K + Ak) Yp - O 


( 54 ) 


where Y 


“ a ^ i * X \ lt 8hould * -entioned that a nonsingular matrix 

d2ter»in^ifn t ^"^I! ed i? Ptl " a J ly by applying coranon mathematical tools for the 
determination of the linear independence of a given number of vectors as 

£ in r€fe ' ence 8* These methods are also^ppU^aWe 
Practical W aDoii e ” u “ bec of constraints higher than the rank of matrix Y 0 . 
reference g™ structu ral dynamics problems are presented in 

«— »• — thi.isjs.1 

ProceduL C ° U ii%^ e contr °^ s to the main structure using the above described 
liHI? doing so, the actual nonlinear stiffnesses C e are replaced by 
linear stiffnesses chosen to be an average representative of ‘the nonUnear 7 

< 2) Caiculate the normal mode characteristics of this linearly coupled 
iJl Setei? Cat * the Unsteady aerodynamic forces based on this set of nor- 

culatims^tteV^M °J, hin9 ® stiffn «ss variations or nonlinear flutter cal- 
be used?** of concepts III and II described subsequently may 


Combined Application of Concepts I, II, and III 
examination of the possibilities offered by the three concepts 
ciai. Four possible variations can be outlined as follows: * 

Combinati on of Concept ill and Concept II : 

hln* SSt.”*’ tMlm ’ lnt ° aCC °“ nt and U ’ htl r 

system^ Caloulate the normal mode characteristics of the linearly coupled 
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Combination of Conc ept Hi and Concept T ; 


trol h «.*„S?L 1 r 1]r ' l H* CO “ PU " 9 *■» — 

controls. rcsultins in s configuration „ith rigidly locked 


* sraiarrsi 

r* j . » _ 


Combination of Co ncept II and Concept t « 


configuration. the ,ICCt!,ft Stcuc "“ e “ lo > controls r.novsd as a basic 

ioc JVSfiiS UppSSS ^ lti ri t ;r° n ,r‘ th “fwiy 

adding modal mass coupling matrices t-o th** 11 * . Thls can ** achieved by 

configuration similar to those defined in equations iT^ ^ ° f ^ baSiC 

cooplSr^tr^Lf^ia'”^ 1 ^^^ co.«ici„ts cl tbs bass 


^rs “ *Rr T ^ A R*Rs 


(56) 


where 


*Rr T * fxr* u yr , u 2r , nxrr T) yi , t) zr ] 


(57) 


and rotationai «•*-- 

system (see fig. 10) . if the center ucture in relation to the XYZ axis 
coupling point, (x s j y_, Za ) - ( 8 o °f the contro1 lies outside the 
written in the form *' °' z) * th ® inertia matrix Aa r can be 


®R 


*** 


| m R®Z 
0 


*"R*z 

0 

®R®X 


®R 


" m R®x 

0 


-m p s 2 

0 


m R 8 z 

0 

“ m R*x 


Irx + ®rb z 0 

0 *Ry + »r(®£ + s|) 

_m R B x 8 z 0 


0 

m R®X 

0 

“ n 'R 8 z®x 

0 


( 58 ) 


*Rz + ®R®x 


where 





•g mass of the control surface , 

I Rx t i Ry , i Rz mass moments of inertia of the control surface in 

y relation to its center of gravity 

(3) Take into account the control degrees of freedom according to Concept I 
by adding a separate set of control normal modes with the main structure at rest. 

Conhination of Concept II and Concept III : 

(1) Test the aircraft structure with rigid control dummies in locked con- 
dition as a basic configuration. The rigid dummies are used to determine a bet- 
ter basic set of normal mode shapes representing the dynamic deformations of 
the coupled system than can be determined in the test configuration with 
removed controls. This procedure can best be described as convergence accel- 
eration by means of interface loading. 

(2) Establish analytically a second configuration with the dummy controls ' 

removed. This can be achieved in accordance with Concept II by subtracting a 

modal mass coupling matrix Am as defined in equation (56) from the equations 
of motion of the basic configuration. 

(3) Apply Concept III coupling the elastic controls to the main structure. 


COMPARATIVE CONSIDERATIONS 

The concepts presented offer a number 3f possibilities to incorporate the 
control systems of light airplanes, which in general are affected by strong con- 
centrated nonlinearities, into the flutter analysis. Special emphasis is placed 
on the mathematical modeling of the elastcmechanical system based on GVT. It 
is obvious that a final evaluation of the applicability and accuracy of the dif- 
ferent concepts is rather difficult because, up to the present time, only 
Concept I has been applied to some extent to real airplane structures. Only 
little experience with the other concepts is at hand. Thus, Concept II has 
recently been employed in the course of the flutter clearance process of the 
soaring airplane ASW-15. Flutter calculations based on this concept predicted 
tail flutter at about 200 km/hr. That result was verified by flight flutter 
tests, where the airplane showed nonlinear flutter in a speed range from 
175 to 220 km/hr, starting with comparably small amplitudes at 175 km/hr and 
increasing to very high amplitudes far beyond the regular rudder stroke at 
higher speeds. This behavior in concurrence with substantial alterations of 
the flutter modes is symptomatic of highly nonlinear flutter cases. A more 
detailed consideration of this special problem exceeds the subject of this 
paper and should be reserved for further investigations. 

It is also worth mentioning that the ground vibration test carried out in 
accordance with this concept toc>. far less test time than a normal test on the 
unchanged structure (reduction of about 80%). 
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The first comparative investigation of the Concepts I, II, and III has 
been the special concern of reference 10, where results are reported for a 
simple plate-type wing-aileron model with largely linear elastodynamical prop- 
erties. Although this model cannot be considered representative in all respects 
of the elastodynamical behavior of real airplanes, it seems to be opportune to 
use the results of this investigation together with the present experience with 
the Concepts I and II as a basis for a preliminary assessment concerning the 
advantages and the weak points of different methods. For this purpose a selected 
number of criteria is used taking into consideration several requirements such 
as 


(1) Test effort required 

(2) Numerical effort required 

(3) General applicability 

(4) Physical consistency 

Table 1 shows in a condensed form how the criteria are met by the several 
concepts. 


CONCLUDING REMARKS 

It has been known for many years that the flutter clearance of light air- 
planes can be highly afflicted by uncertainties stemming from strong localized 
nonlinearities in the control mechanisms. It is shown that the establishment 
of more reliable and accurate mathematical models for the flutter analysis 
requires modified ground vibration test procedures combined with suitably 
adapted modal synthesis approaches. Three basic concepts with several varia- 
tions have beer described in detail. They offer a diverse choice of tools for 
carrying out both approximately linearized and nonlinear flutter investigations. 

A comparative consideration has been made as to the capacity as well as 
the drawbacks of the different concepts. However, because of lack of practical 
experience with Concepts II and III, it is not possible at present to make a 
conclusive evaluation. 
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Figure 1.- Schematical sketch of the control system of a light airplane. 



(a) Force deflection diagram. 



0. deg 


DAMPING 

DECREMENT, 


C e '(3> 


(b) Hinge stiffness and damping versus hinge angle. 


Figure 


2.- Force deflection diagram and stiffness and 
for the rudder system of a soaring airplane. 


damping 
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AILERON HINGE ANGLE, 
fc deg 


(a) Force deflection diagram. 



AILERON HINGE ANGLE. 

3. deg 

(b) Hinge stiffness and damping versus hinge angle. 

Figure 4.- Force deflection diagram and stiffness and damping for 
the aileron system of a soaring airplane. Symmetrical case. 




Cj = 4.7 N-m 
C ? = 16 N-m 


Figure 5.- Measured flutter boundary of a nonlinear model. 



NORMAL MODE WITH 
CONTROL SURFACE LOCKED 


NORMAL MODE WITH LIFTING 
SURFACE AT REST AND 
CONTROL SURFACE FREE 


ARBITRARY DEFLECTION 


Figure 6.- Modal superposition according to concept I 
Controls without tabs. 




Figure 7.- Lifting surface with control surface and tab. 


0 NORMA' MODE WITHCS AND 
Ar T LOCKED 






® B NORMAL MODE WITH CS FREE 
T LOCKED. ANDES AT REST 


fl>, 


Cv NORMAL MODE WITH CS LOCKED. 
T FREE. AND LS AT REST 



ARBITRARY DEFLECTION 

LS: LIFTING SUR : ACE 
CS: CONTROL SURFACE 
T: TAB 


Figure 8.- Modal superposition according to concept I, 
Controls with tabs. 


465 



4 

^ 79-27079 

ADVANCED COMPOSITES in SAILPLANE STRUCTURES: i 

APPLICATION AND MECHANICAL PROPERTIES 
Dieter Muser 

Research Center Stuttgart 
Deutsche Forschungs- und versuchsanstalt 
fur Luft- und Raumfahrt e.V. 


SUMMARY 


Advanced Composites rangers to 18% in compar- 
ers in an epoxy matrix. Weigh 9 laminates will be produced by hand-layup 

Ison with glass fib.. structures. ™ done with th... mete- 

technlques and all material calculation of strength and stiffness as 

rials. The" value. “ y lals to get a weight-optimum construction. 

^irfor^ri-:Uu» construction are mentioned. 

technical history 


The first fiber-reinforced 9 J^ er ^ e than 4000 glides with glass-fiber- 
made its maiden flight in 195 • ' u t the world. Increasing the wing 

reinforced structures are but structural demands increased 

loading permitted increases in maximum speeo, 

the weight also* 

A large span enabled the 5°" St '^us 8 2-° 49) ^n^ sinking speeds of 0.50 m/s 
ratios of about 50 (ASW 17: 48 * 5 ’ *; Uz * w ing spans with more than 

( 1.«4 ft/s). But it was not poms ^/possible when carbon 
22 meters without a very soft w ' * section of the Akaflieg Braunschweig SB 10 

fibers were used in the center "infection or glider ha8 the 

in 1972 (fig. D. With o iTs 0.35 ft/s). But the 

best glide ratio of 53 and a sinking speed o ^ so this mate rial was used 

price of carbon fibers was very high at thi t fs _ 29 in 1975 . Xo realize 
only in another prototype, the Akaflieg 9 absolutely necessary to 

the old dream to vary the span during night, it wa ^ ^ ^ ^ of the 

use carbon fibers in the outer .^ml/Traunschweig built the first all-carbon 

.“**■ - u “ d ln ■ ,ar, “ 

ume in different types of commercial gn 
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